Giant Orbital Hall Effect in Transition Metals: 
Origin of Large Spin and Anomalous Hall Effects 
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In transition metals and their compounds, the orbital degrees of freedom gives rise to an orbital 
current, in addition to the ordinary spin and charge currents. We reveal that considerably large 
spin and anomalous Hall effects (SHE and AHE) observed in transition metals originate from an 
orbital Hall effect (OHE). To elucidate the origin of these novel Hall effects, a simple periodic s-d 
hybridization model is proposed as a generic model. The giant positive OHE originates from the 
orbital Aharonov-Bohm phase factor, and induces spin Hall conductivity that is proportional to the 
spin-orbit polarization at the Fermi level, which is positive (negative) in metals with more than (less 
than) half-filling. 

PACS numbers: 72.25.Ba,72.10.Bg,72.80.Ga 
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The Hall effect, first discovered at the end of 19th cen- 
tury, has revealed the profound nature of electron trans- 
port in metals and semiconductors via the anomalous 
Hall effect (AHE) and (fractional) quantum Hall effects. 
It has recently been recognized that conventional semi- 
conductors and metals exhibit a spin Hall effect (SHE), 
which is the phenomenon where an electric field induces 
a spin current (a flow of spin angular momentum s z ) in 
a transverse direction [1-5]. Recently, a theory of the 
intrinsic Hall effect proposed by Karplus and Luttinger 
[6], which occurs in multiband systems and is indepen- 
dent of impurity scattering, has been intensively devel- 
oped [7, 8]. In particular, a quantum SHE has also been 
predicted and experimentally confirmed [9, 10]. 

The spin Hall conductivity (SHC) observed in transi- 
tion metals has given rise to further issues regarding the 
origin of the SHE, since the SHC observed in Pt exceeds 
200 he" 1 ■ fi — 1 cm _1 , which is approximately 10 4 times 
larger than that of n-type semiconductors [5], and the 
SHCs in Nb and Mo are negative [11]. The large SHE 
and the sign change of the SHC in transition metals has 
attracted much interest, and many theoretical studies of 
the SHE have so far been conducted based on realistic 
multiband models for Ru-oxide [8] and various Ad and 
5d metals [12], including Au, W [13], and Pt [14, 15]. 
The calculated results for the SHC semi-quantitatively 
agree with the observed results. The mechanism for the 
SHE has been explained in such a way that spin-orbit 
interactions (SOI) and the phase of hopping integrals of 
electrons give rise to the Aharonov-Bohm (AB) effect, 
and therefore the conduction electrons are subject to an 
effective spin-dependent magnetic field. 

Since the transition metals have orbital degrees of 
freedom in addition to the spin and charge degrees of 
freedom, flow of the atomic orbital angular momentum 
(l z ), that is, an orbital current, may be realized in a 
nonequilibrium state. In fact, several authors have pre- 



dicted the emergence of a large orbital Hall effect (OHE) 
[8, 12, 16], which is a phenomenon where an electric field 
induces a flow of p- and d-orbital angular momentum in 
a transverse direction. In particular, the predicted or- 
bital Hall conductivity (OHC) in transition metals and 
oxides [8, 12] is considerably larger than the SHC. Figure 
1 shows the OHC and SHC calculated for transition met- 
als using the Naval Research Laboratory tight-binding 
(NRL-TB) model [17]. In each metal, the magnitude of 
the OHC exceeds the SHC, even in topological insulators 
(e.g., ~ e/27r in graphene [9] and HgTe [10]). Interest- 
ingly, while the SHC consistently changes its sign with 
the electron number n, as with several recent experiments 
[5, 11], the obtained OHC is almost independent of the 
SOI and is always positive. These prominent and univer- 
sal features of orbital dynamics in metals, independent 
of crystal and multiband structures, have not been rec- 
ognized until recently. 

In spite of these remarkable features of the SHC and 
OHC given in previous works [8, 12], no physical origin 
of the "giant OHE" nor "hidden relationship" between 
OHE and SHE have been presented. We will show be- 
low that the key phenomenon is the orbital Hall current, 
which originates from the "orbital AB phase factor" that 
reflects the phase factor of the c?-orbital wavefunction. 
Then, the SHC is approximately given by the product 
of OHC and the spin-orbit polarization due to the SOI. 
Even the AHE is understandable in the same concept. 

In this Letter, we discuss these intrinsic Hall effects 
in a unified way by proposing a simple s-d hybridization 
model as a generic model, and explain why the OHC is 
positive and much larger than e/2ir in each transition 
metal. We stress that the large SHE in transition met- 
als originates from the OHE in the presence of atomic 
SOI, not from the Dirac point monopole as in semicon- 
ductors. The derived SHC is approximately proportional 
to the spin-orbit polarization, which is positive (nega- 



2 



tive) in metals with more than (less than) half-filling, 
which is consistent with recent experimental observations 
[4, 5, 11]. It is noted that the present OHE is differ- 
ent from the Hall effect of the angular momentum of s- 
electrons r x p discussed in Ref. [18]. 
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FIG. 1: (Color online) SHC and OHC obtained in ref. [12]. 
The body-centered cubic structure for n = 5, 6, the hexagonal 
closed packed structure for n = 7, 8, and the face-centered 
cubic structure for n = 9 ~ 11. n<i ~ n — 1 in each metal. 
The SHC in Pt is approximately 1 in the unit of e/2ira ~ 10 3 
f2 _1 cm _1 , where a — 0.39 nm is the lattice constant in Pt. 
R = (I z s z )fs is the spin-orbit polarization for 5d metals. 

In fid transition metals (n — 3 ~ 5 is the main quan- 
tum number), the electronic states near the Fermi level 
(/i) are constructed by fi<i-orbitals, which hybridize with 
free-electron-like (n + l)s- and (n + l)p-bands, and form 
a narrow band with a bandwidth in the order of 1 eV 
[17]. The band structures of the transition metals are 
well characterized by the lattice structures. However, the 
OHC seems to be independent of the atomic species of the 
transition metals as shown in Fig. 1. The result suggests 
that the details of the band structure may be irrelevant 
to the occurrence of the large OHC in transition metals. 
Therefore, two-dimensional (2D) s-d hybridization model 
with the orbital degree of freedom was adopted to study 
the universal nature of the OHE in transition metals. 
The s-d hybridization model may also be the simplest 
version of the muffin-tin potential approximation, where 
localized d-electrons can move only via s-d hybridization. 

The Hamiltonian for the 2D s-d hybridization model 
is given by H = J2k <y ai a H a k(T , where 



Hn = 



kL 



k-L 



E d 









(1) 



and *a k(T = (c k(T , d kLa , d k -Lcr); c k(T and d kM a are annihi- 
lation operators for the s- and d-electrons, respectively. 
a = ±1 (or |, |) is the spin index, and M = ±L (L is 
a positive integer) represents the angular momentum of 
the d electron. e k and Ed represent the unhybridized en- 
ergies of the s- and d-electrons, respectively, and VkM is 



the s-d mixing potential. In transition metals, the OHE 
and SHE are mainly caused by interorbit transitions be- 
tween d yz and d zx (l z = ±1) orbitals, and d xy and d x 2_ y 2 
(l z = ±2) orbitals [12, 14]; the former (latter) contribu- 
tions can be obtained by letting L = 1 (L = 2) in the 
present 2D model. 

To elucidate the universal properties of the OHC and 
SHC that are independent of the detailed crystal struc- 
ture, we assume that e k = k 2 /2m (plane wave). As is 
well known, e akr = J2n i n J n (kr)e in ^"-^\ where J n 
is the Bessel function, tp^ — tan -1 (k x /k y ) and ip r = 
tan -1 (x I 'y) . Since the wavefunction of the d-electron is 
£m(i") = f(r)e lMVr , the s-d mixing potential is given by 



y Qe iM Vk in the extended 
(Note that V kM -> ^ 



kM 



VkM = JV kr )*ff a/(r)dr = 
Brillouin zone scheme [19-21] 
under the particle-hole transformation.) Here, the k- 
dependence of Vq is neglected. We will show that the 
phase factor of Vum plays an essential role in the OHE, 
and the derived OHC takes a large value irrespective of 
the nonconservation of M = ±L. 

The present model is similar to the periodic Anderson 
model, which has been intensively studied as an effective 
model for /-electron systems [19-21] . The close similarity 
between the periodic Anderson model and the d-p model 
has been indicated in the previous study of the AHE 
[22]. As is well-known, the hybridization band of Eq. (1) 
is given as E± = ±[(e k + E d ) ± ^(e k - E d ) 2 + 2V 2 }, the 
bandstructure of which is shown in Fig. 2(a). In the 
metallic state, the Fermi level fi is located in the upper 
(E£) or lower ) branch, and the relationship (fi — 
£fc F )(Ai - E d ) = 2V 2 holds. The relation N d (0)/N s (0) = 
2Vq /(p — Ed) 2 3> 1 is satisfied in transition metals, where 
N B (0) = m/2-K is the s-electron density of states (DOS) 
per spin in 2D, and N d (Q) is the d-electron DOS at fx. 

According to the linear-response theory, the intrinsic 
Hall conductivity is given by the summation of the Fermi 
surface term (J-term) and the Fermi sea term (//-term) 
[23] . In previous work, we have shown that the /-term is 
dominant in many metals [12, 14, 24]. Using the 3x3 
Green function G k (w) = (ui + fi — Ho)^ 1 , the /-term of 
the OHC at T = is given by [8, 12] 



(2) 



where J y is the y-component of the charge current, which 
is given by —edHo/dk y . In the present model, 
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(3) 



where — e is the charge of the electron. In addition, J® in 
Eq. (2) is the a;-component of the orbital current, which 
is given by J° = { J° , l»}/(-2e) [12], where {l z ) %t3 = 
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FIG. 2: (Color online) (a) Schematic band structure of the s-d 
hybridization model in Eq. (I), (b) Examples of the clockwise 
motions of electrons along the nearest three sites, which give 
the orbital AB phase factor e tLe * with 8i > 0. (c) Localized d- 
orbital state and the Fermi level of conduction electrons under 
E y . (d) A semiclassical explanation for the Hall current (_L E) 
due to the angular momentum conservation. An left-moving 
electron is converted to right-moving by l z = +L orbital. 
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Here, the relations 9Vki/9fc x 



-iL(k y / k^V^L and 



dV^L/dky = iL(k x /k 2 )V k L are used [21]. Thus, the mo- 
mentum derivative of the s-d mixing potential gives rise 
to an anomalous velocity that is perpendicular to k: The 
OHC is proportional to ((J a c ) M ( J° )i, 2 V r k * t ) F s + 0. 

By inserting Eqs. (3) and (4) into Eq. (2), the OHC 
in the present model is simply obtained as 



O z 

w x y 
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mN ^ 
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where N is the number of k-points, g\s.{to) = (lo + /i — 
E£)(uj + fi — and 7 is the damping rate (due to 

impurity potentials). The obtained OHC is finite even if 
SOI is absent [8, 12, 14]. Since lim 7 ^ +0 j/{x 2 + -f 2 ) = 
ttS{x) and jr 8(n - E±) = N d (0), We obtain that 



O: 



eL2 
Am 



N d (0) 



2V 2 



E k ) 2 



47T ' 



(6) 



where we used that the term in the bracket is ~ N s (0) 
when N d (0)/N s (0) > 1. Thus, the OHC takes an ap- 
proximate universal positive value, independent of the 
model parameters. For the 3D s-d hybridization model, 



VkAf is proportional to the spherical harmonics with 
L = 2. In this model, we obtain 0% y ~ (e/2na) ■ ir if 
we put k~F ~ 7r/a, where a is the lattice spacing. This 
result is qualitatively consistent with the results in Fig. 1 
for n < 9. (OHCs in Au and Ag are small since n d = 10.) 
We have verified that the //-term is negligibly small. 

The reason why the giant OHE emerges in transition 
metals based on a 2D tight-binding model with three 
orbitals (s and g?±l) at each site is now discussed. Long- 
range hopping integrals must be considered to reproduce 
a free-electron-like s-electron dispersion [17]. Figure 2(b) 
illustrates two examples of the clockwise motion of elec- 
trons along the nearest three sites [dr, — » s — » s — » di]; as 
discussed in Ref. [12], these are important processes for 
the OHE, because the s-d hopping integrals are much 
larger than the d-d hopping integrals. Therein, the 
electron acquires the phase factor e lLe due to the an- 
gular dependence of the mixing potential in real-space, 
Vi(r) oc e lLVr , where 9 is the angle between the incom- 
ing and outgoing electron (8 — 7r/4 or 7r/2 in this fig- 
ure). This can be interpreted as the "orbital AB phase" 
given by the effective magnetic flux tp ■ (L9/2ir) through 
the area of the triangle, where (f>o = 2nh/e is the flux 
quantum. It is simple to check whether any of the other 
three-site clockwise motions cause the factor e lLe with 
9 > 0. Therefore, the d-electron with l z = ±Z is sub- 
ject to the huge and positive effective magnetic field that 
reaches ~ ±4>o/a 2 . This is the origin of the giant positive 
OHE in the order of e/27r in transition metals. 

We also present a simple semiclassical explanation, in 
that the OHE is a natural consequence of the anisotropy 
in s-d hybridization under the electric field E y . In the 
muffin-tin model, the localized d-electron (virtual bound 
state) moves to the conduction band via the s-d mixing 
potential [25] as shown in Figs. 2 (c) and (d). According 
to the Fermi's golden rule, the d — >s tunneling probabil- 
ity at point B is given by T B ~ J^ g duN d (u)\V \ 2 N s (u), 
where (1b = ^ — e.E y r d is the electrochemical potential at 
B under E y Then, Ta — —^b and ~ eE y r d since 
N s (0) ~ l/i-e^l" 1 and N d (0) ~ |/i - E d \-\ Be- 
cause of the angular momentum conservation, the ve- 
locity of emitted (absorbed) electron at point B (A) 
has rr-component; v% ~ (—)L/r d m. If we assume that 
the tunneling s-electron hybridizes to one of the neigh- 
boring sites and l z = ±L is quasi-conserved, the cur- 
rent of the successive tunneling electron will be ~ 
inlr^sl • asin^o, where n is the electron density and 
sin(f = \v%\/vf ~ O(l). Therefore, the estimated trans- 
verse ri-orbital current density is 



J., 



^2 M if ~ eL 2 E y na/mv F . 



(7) 



M=±L 



Since a ~ 7r/c F 1 and n 



a 2 , 0% y = j° ' /E y is in the 



order of +eL 2 . Thus, Eq. (6) is reproduced (aside from 
a numerical factor) by this semiclassical consideration 
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We note that the partial wave of the l z = L chan- 
nel, tPl{y) cx JL(kr)e tLVr , has a small overlap inte- 
gral between the nearest sites, / ^(r)^!/ (r + ax)dr, for 
V = —L, due to the phase factor in iPl(y)- For this rea- 
son, l z is quasi-conserved when the tunneling s-electron 
hybridizes to the rf-orbital at a neighboring site. Mathe- 
matically, the anomalous velocity is given by taking the 
gradient of the phase factor in Vi( r )j see c 1- (4)- 

Next, we discuss the SHE in the presence of the atomic 
SOI; \ ■ ^ (A > 0). Since (M\l v \M') = for v = 
x, y in the present 2D model, the atomic SOI for the o- 
spin is given by (#£)ij = (Acr/2)L(<5;, 2 - $i,3)6i,j- Also, 
only the z-component of SOI is significant for the SHE in 
real transition metals [12, 14]. Using the Green function 
G kCT (w) = (uj + h-H - the SHC is given by 

1 .V^Tr 



2irN 



k,cr 



J, 



c Gka(n)ifG kff (-n) (8) 



If the A-dependence of eigenenergies is neglected, which 
corrects the SHC of order 0(A 3 ) to, Eq. (8) becomes 

a z xy « 2R/L 2 



■01 



(9) 



where R = (l z s z )ps represents the spin-orbit polarization 
ratio due to the SOI at the Fermi level, which is given 
by R = LX/(p — Ej) in the present model up to 0(A). 
Thus, the SHC is positive (negative) when /i is located 
in the upper branch E£ (lower branch E^). 

It is natural to expect that the relationship in 
Eq. (9) holds in real transition metals, where 
the spin-orbit polarization ratio is defined as R = 
E m / FS (^S2)k,mrfSk, m / E m / FS rfS k,m in real systems, 
where m is the band index. To verify this expectation, 
R is shown for 5d metals given by the NRL-TB model in 
Fig. 1: The obtained R is positive (negative) in metals 
with more than (less than) half-filling, which is consistent 
with Hund's rule. The qualitative similarity between the 
SHC and R in Fig. 1 indicates that the spin current j s 
is induced by the orbital current j in proportion to R, 
and therefore the relationship in Eq. (9) holds approx- 
imately for various metals. (In fact, Xl z s z provides the 
dominant contribution to the SHE [12].) As a result, 
the present analysis based on a simple s-d hybridization 
model captures the overall behavior of the OHE and SHE 
in transition metals. 

In the low resistivity regime, the intrinsic SHC is given 
by integrating the k-space Berry curvature of Bloch wave- 
function (Berry curvature term) [1, 2, 12, 14, 15]. In 
fact, previous studies based on the tight- binding mod- 
els [12, 14] and the band calculation [15] had succeeded 
in reproducing experimental SHC's in several transition 
metals with low resistivity, both in magnitude and sign 
[5, 11]. The present study has shown that the large Berry 
curvature in transition metals, the origin of which had 
been unclear in previous studies, originates from the d- 
orbital angular momentum. We have revealed the exis- 
tence of the real-space orbital Berry phase (— AB phase) , 



which causes not only the giant positive OHE without 
using the SOI, but also the SHE and AHE if R ^ 0. By 
virtue of this scheme, hidden relationships between the 
OHE and SHE have been derived. Although the OHE is 
indirectly observed via SHE and AHE, it is interesting 
to detect the OHE directly: We propose that a meso- 
scopic "H-shape" circuit will be useful, which was origi- 
nally used to measure the SHC in scmimctals [26]. 

To summarize, we have revealed that the giant positive 
OHE in transition metals originates from the orbital AB 
phase due to the d-orbital angular momentum, without 
necessity of any special bandstructure (e.g., Dirac point 
monopole at n). We have shown that the OHE is the es- 
sential phenomenon, and it induces the large SHE (AHE) 
in paramagnetic (ferromagnetic) metals in the presence 
of SOI. The sign of the SHC is equal to that of the spin- 
orbit polarization (Hund's rule), which is consistent with 
recent experimental observations [4, 5, 11]. An intuitive 
explanation for the intrinsic Hall effect in real-space is 
presented in Figs. 2 (c) and (d). 
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been supported by Grants-in-Aids for Scientific Research 
from MEXT, Japan. 
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